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Abstract 

Hybrid meson energies are calculated in the static quark limit with the Dy- 
namical Quark Model (DQM). In the DQM, transverse gluons are represented 
as effective constituents with a dynamically generated mass. Hybrid masses 
are determined within the Tamm-Dancoff approximation for the resulting rel- 
ativistic Salpeter equation. Although the general features of the adiabatic 
potential surfaces correspond with lattice data, the results disagree on level 
orderings. Similar problems appear to exist in all constituent glue models of 
hybrids. We conclude that constituent gluons do not accurately represent soft 
gluonic degrees of freedom. The steps necessary to correct this deficiency are 
discussed. 
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I. INTRODUCTION 



A decade of experimental signals [1IJ for QCD hybrids (in particular with J PC = 1 h ) has 
culminated in the claimed observation of three such states at Brookhaven [0 in the last year. 
The question of the nature of QCD hybrids has thus become topical. Furthermore, lattice 
gauge calculations are now at the point of accurately determining light hybrid masses. In 
view of these developments, it is of interest to compare models of strong (low energy) QCD 
with lattice data to determine their viability and to explicate and guide current experimental 
efforts. 

It is often stated that a hybrid is a hadron consisting of valence quarks and glue. However, 
one must specify what is meant by the notion of valence glue for this statement to be useful. 
There are two broad ideas in this regard: it is some sort of string or flux tube ||,|J or it is 
an effective constituent confined by a bag or potential ]7|-|§ . 

As an example of the importance of choosing correct degrees of freedom, we mention 
the simple problem of determining the number of components of a constituent gluon. It 
has been suggested that a massive constituent gluon should be transverse so as to maintain 
consistency with Yang's theorem ||. However it was noted that this is inconsistent with 
the requirements of Lorentz invariance. Thus, for example, J = 1 glueballs are expected to 
exist and lattice calculations indicate that they are quite heavy (roughly 3 GeV) [Tt|. Such 
a state may not be constructed from two transverse constituent gluons (Yang's theorem) 
and therefore may be expected to have a mass of roughly 3m g ~ 3 GeV. However massive 
vector gluons have no such constraint and one therefore expects them to have a mass of 
approximately 2m g ~ 2 GeV. 

The nature of the appropriate effective degrees of freedom for glue can only be deter- 
mined by a long process of calculation and comparison with experimental and lattice data. 
There are, however, a few indications that low energy glue is string-like. Perhaps the most 
compelling of these are lattice calculations of energy, action, or field densities between static 



color sources which are reminiscent of flux tubes |11|. An intriguing clue is also provided 
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by the spin splittings of heavy quarkonia. It is known that an effective interaction free of 
long range exchange spin-orbit coupling is needed to reproduce the mass splitting of the 3 Pj 
heavy quarkonium multiplets. However, an analysis of QCD in the heavy quark limit con- 
vincingly demonstrates that obtaining such an effective potential requires that low energy 
glue must be string-like [|T^ . 

Alternatively, pointlike models of low energy glue have a long history, originating with 
MIT bag model calculations of Barnes f| and others ||. Horn and Mandula @ were the 
first to consider a potential constituent glue model of hybrids. Their hybrids consisted of 
constituent quarks and pointlike, massless, spinless, and colorless glue in a nonrelativistic 
potential model. The confining potential was taken to be linear with a string tension given 
by the ratio of color Casimir operators: b qg = 9/8b qq -. The authors noted that the two body 
qq potential is anti-confining in the color octet channel and has a repulsive Coulomb spike at 
short distances. They argued that this is unphysical and hence choose to neglect this term 
in the interaction. It is clear that a great many simplifying assumptions have gone into the 
construction of this model. It is our purpose to compare a more sophisticated version of the 
model to lattice data to learn something about these assumptions. 

In the following we employ a model field theoretic Hamiltonian of QCD. The model incor- 
porates linear confinement at low energy and evolves into perturbative QCD at high energy. 
A nontrivial vacuum is used to generate constituent quark and gluon masses. The eigenvalue 
equation is derived for a qqg system where the quarks are static. The resulting adiabatic 
potential surfaces are then compared to recent lattice results. We conclude that the simple 
picture of glue as a pointlike constituent particle reproduces the general behaviour of the 
lattice results but fails to yield the correct level orderings. Furthermore, other models which 



regard the gluonic degrees of freedom as pointlike (eg., fl^fl , [[/J) do not contain sufficient 
degrees of freedom to generate all of the adiabatic surfaces. Thus constituent glue models 
appear to fail to describe hybrids. This stands in contrast to string or bag-like models which, 
although disagreeing on details, capture the rough structure of the lattice data. 
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II. A CONSTITUENT GLUE MODEL OF STATIC HYBRIDS 



A. The Dynamical Quark Model 

The starting point for our description of hybrids is the following model Hamiltonian: 

H = J rfx^(x) [-ia ■ V + (3m] ^(x) + \ J [|E A (x)| 2 + \B A (x)f 

+ i/dxdy/(x)F(|x-y|)/(y) (1) 

where the color charge density is p A (x) = ^^(x)T j4 '0(x) — f ABC A B ( y x) ■ E c '(x) and the 
potential is given by 

"M = T-]Pr( 1 -<>- w ) w 

and N c — 3. The quark mass appearing in this Hamiltonian is the current mass. To be 
phenomenologically successful, constituent quark masses must be generated in some way. 
This may be achieved by employing a BCS vacuum Ansatz; the gap equation which follows 
from minimizing the vacuum energy density, (Q\H\Q), (where |0) represents the BCS trial 
vacuum) gives rise to a low energy constituent quark mass of roughly 200 MeV |13[ . A similar 
calculation in the glue sector yields a gluon dispersion relation which is well-approximated 
by 

oo(k) = ^k 2 + m 2 g e- k / K (3) 

with m g = 800 MeV and k = 6.5 GeV. One sees that a constituent gluon mass of approxi- 
mately 800 MeV has been generated [Q. Hadrons are then constructed on top of the BCS 
vacuum \Q) by employing a basis truncation (typically Tamm Dancoff or Random Phase) 
and solving the resulting Bethe-Salpeter equation. This approach has been used to derive 
the low lying spectrum of glueballs and agrees remarkably well with lattice data ||14|| . It 



should be noted that the dynamical gluons are transverse so that Yang's theorem holds and 
the difficulty mentioned in the Introduction does not arise. 
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In the following, the parameters of the model are fixed to the qq potential derived in the 
lattice calculation of Juge, Kuti, and Morningstar |L5] (see Fig. 3 below). The fit yields 
a s = 0.29 and b = 0.24 GeV 2 . The final parameter, Auv serves as an ultraviolet cutoff on 
the linear confinement potential. Its value was determined in Ref. JTJJ by fitting the gluon 
condensate, and will be set at 4 GeV in the following. We shall consider the static quark 
limit in the remainder of this work so that quark masses (and the quark sector of the BCS 
vacuum) become irrelevant. 

The model presented here may be considered as a simplified version of the Coulomb 
gauge QCD Hamiltonian where the effects of nonabelian gauge couplings have been modeled 
by the linear confinement term. Furthermore, second order transverse gluon exchange is 
suppressed by the heavy quark masses. In principle, this approach allows the elimination 
of the ultraviolet scale Ajjy v i & renormalization. A method for achieving this which is 
appropriate for nonperturbative Hamiltonian-based calculations is described in Ref. ||13|| . 



B. Static Hybrids 

In the Tamm Dancoff approximation, hybrids are constructed as qqg excitations of the 
BCS vacuum. For the heavy hybrids considered here, the (anti) quarks serve as static color 
sources (sinks) and the gluons are constituent particles as described above. We choose to 
work in the "diatomic molecule" basis because this facilitates comparison with the lattice 



results of Ref. |15[]. Thus the hybrid state may be written as 

\Rn g ;j g A0 = J dk^ ng3g (fc)2?ft {R)V^ (k) ^ ^^xi v ^^/.r^Iax I fi > • ( 4 ) 

Wigner rotation functions are written as T> in this expression, R is the distance between 
the qq pair, A is the gluon helicity, and (p is the radial hybrid wavefunction in momentum 
space. The gluon polarization wavefunction, X\\> is given by x\y — fi\x/V2 an d X\\> = 
\5\\r/^/2. The two cases, £ =1, -1, represent transverse magnetic and transverse electric 
hybrids respectively. Finally, A is the projection of the gluonic angular momentum onto 
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the qq axis, j g is the total angular momentum of the gluon, and n g labels the radial basis 
state. We note that employing helicity basis gluon creation operators makes this expression 
significantly more compact than the canonical basis. 

The Salpeter equation which follows from this ansatz and the Hamiltonian in Eqn. (1) 
may be obtained from the following matrix element: 



(R'n' g ,f g A>Z'\H\Rn g]Jg AO = J dk ip* n , fg (k)p ngjg (k)-[u(k) + — - 
+ 1 ( [dqdUhdk)^ gjg (k)V(k-q) + ^ 



a + (q-k) 2 ) 



uj(k)ui(q) 

dqdk V l, rg {q) Vngjg {k) V(k - q) (e<#<*-«> + e^^) ^=== X 

'u(k)ui(q) 



+ (qV(r) + -y(o)) K g f g n ajg Sf gj M^ (5) 

An overall 5(R' — R) is understood and Nri g j' g n g j s = J dkk 2 ip* n ,^{k)ip ng j g {k) is the wavefunc- 
tion normalization factor. The two extra Wigner rotation matrices arise from converting 
the Cartesian basis implicit in Eqn. (1) and in the expression for the color current to the 
helicity basis, a^A = ^H{kX)dkXA = ^mxi^y ^ — 0) e c( m ) a fcAA, where 6q and e# are canonical 
and helicity polarization vectors respectively. 

The first term in this expression is the gluon kinetic energy, the second is the gluon self 
energy, the third is the gluon potential, and the fourth is the qq potential and self energies 
for static quarks in a color octet. The presence of the gluon and quark self energies assures 
the infrared finiteness of the Salpeter equation. This appears to be a general feature of color 
singlet states in our approach [[b| . 

The last task is to identify the diatomic quantum numbers used to label the states. In 
the Jacob- Wick convention the action of parity on glue is given by 

P^XA p] = ^e- 2 ^a_^ A (6) 

where <p is the azimuthal angle of k and rj^ — — 1 is the intrinsic gluon parity. Thus the 
hybrid states given in Eqn. (4) are eigenstates of gluonic parity with 



P g \Rn g ; 3g AO = ^(-)^ +1 \Rn g] j g A£). (7) 

A reflection of the glue degrees of freedom through a plane containing the qq axis leaves the 
Hamiltonian invariant and when acting on the states it takes A — > —A. For |A| > one 
thus has doubly degenerate states. We call this operation y-parity and perform it by setting 
R — > Rz and taking y — > —y which may be achieved by a parity operation followed by a 
rotation through 7r about the y-axis. The action of Y on a single transverse gluon state, 
\k\A) 

= a l\A\fy i s therefore given by, 

Y\k\A) = e- i7rJ yP\k\A) = e 2iA *|£;', -A, A), (8) 

where k! = (k x , —k y , k z ). For the states defined in Eq. (|J) one has 

Y\Rn g ;j g AO = ^(-) A+1 \Rn g ;j g - AC), (9) 

where the relations i{k') = (— ) 2j+M+M 'e~ 2v '^X>i ,(A;) and ,(z) = 5^ were used. 
For A / the Y-diagonal states may thus be written as 

\Rn g -J g \A\^ VY ) = ^(\Rn g - 3g \A\i) + 7 lY \Rn g - 3g - |A|£)J, (10) 

where r] Y = ±1 and, 

Y\Rn g ;j 9 \A\Z;r]y) = ^ VY (-) A+1 \Rn g ;j g \A\^r] Y ). (11) 

For A = we simply have, 

Y\Rn g ;j g 0O = -^\Rn g -j g ^). (12) 

To allow for easy comparison with conventions used elsewhere, we summarize the trans- 
formation properties of qqg states under total parity, P = P q P g , and charge conjugation 
C = C q C g , including the quark degrees of freedom. 

P\Rn g ;j g \A\^rj Y ) = £« 9 (-)^ +1 | - Rn g ;j g \A\^ VY ) (13) 



C\Rn g ;j g \A\Z; VY ) = r/^(-) 5 ^ +1 | - Rn g ; j a |A|£; rj Y ) (14) 
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where r]? q = r]^ — rQ = -1. The states introduced in Eq. (|T0D are therefore also eigenstates 
of combined PC, with 

PC\Rn g ;j g \A\^rjy) = ^^{-y^Rn^j^^y) (15) 

After dividing out the quark portion of PC we are left with 

(PC) g \Rn g ;j g \A\Z;r]y) = {-)^ +1 \Rn g] j s |A|£; Vy ) (16) 

C. Results 

The Salpeter equation is solved by expanding the radial wavefunction in a complete 
basis and by diagonalizing the resulting Hamiltonian matrix. The evaluation of the matrix 
elements is greatly facilitated by performing the angular integrals analytically. Furthermore, 
we found it expedient to do the remaining numerical integrals with the potential in position 
space since the integrand is less oscillatory for large argument in this case. The plane 
waves and potential were expanded in a double series of Wigner functions yielding a total 
of ten Wigner functions. The angular integral then evaluates to a sum over a product of six 
Clebsch- Gordon coefficients. 

As discussed above eigenstates of the Hamiltonian may be labelled with the projection 
of the angular momentum onto the qq axis, the product of gluonic parity with charge con- 
jugation, (PC) g and the Y-parity eigenvalue. For A / the two Y-parity eigenstates are 
degenerate. States are therefore denoted by Aj PC ^ g where A = 0, 1, 2 are denoted by S, II, A; 
(PC) g = TjgVg ^(~y 9+1 = g or u for even or odd parity respectively; and Y = C,f]gr]y(—) A+1 
= ±. The total gluonic angular momentum, j g is not a good quantum number. We have 
found however, that the sum over radial basis states and gluon angular momentum satu- 
rates quickly. For example, Fig 1. shows the (approximately) exact wavefunction for E+ 
with j g = 1 and a variational single Gaussian orbital. The next contribution to the E+ 
eigenstate has j g = 3, this is shown as the dotted line in the figure. Evidently the eigenstate 



S 



is dominated by the lowest gluonic angular momentum component and the single Gaussian 
approximation is quite accurate. This remains true for all R studied here, with the j g = 3 
component rising to 13% of the wavefunction for R = 10 GeV -1 . These observations support 
traditional approximations made for hybrids: the use of simple Gaussian wavefunctions [19| 
and the truncation of the adiabatic Schrodinger equation to lowest j g 

FIGURES 




Fig. 1. wavefunctions. Exact result (dashed line), single Gaussian approximation (solid line), 

j = 3 component (dotted line). 



The results for the gluon spectrum are presented as a function of the static quark sep- 
aration in Figs. 2 and 3. They are plotted in terms of the lattice scale, Rq = 2.32 GeV -1 
|i~5l and the potentials have all been normalized by subtracting an overall constant given 
by V q g(2R Q ). In Figs. 2a-c we compare the recent lattice results from Ref. |15| with the 
predictions of the flux tube model Q. The flux tube model was motivated by the strong 
coupling limit of the QCD lattice Hamiltonian. It is based upon a nonrelativistic, small 
oscillation approximation to motion of the colored links in a topological sector where there 
are no overlapping (color representations of higher dimension) or disconnected links. Under 
parity (charge) conjugation the spatial (color) orientation of lattice links is reversed, thus the 
nonrelativistic "beads" of the flux tube model are assumed to flip orientation with respect 
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to the position of the quarks under parity and to have positive intrinsic charge conjugation 
parity. 

The lowest solid line in Fig. 2a is a fit to the qq potential given by E(R) = —4as/(3R) + 
bR + const which corresponds to the ground state lattice potential. The other solid lines 
in Fig. 2 show the flux tube potential as given by 

E (R) = bR + ^-(l-e-f bl/2R ) (17) 
R 

with / ~ 1 and N = J2 m =i m ( n m+ + n m -). The latter represents the total number of right- 
handed (n m+ ) and left-handed (n m _) transverse phonon modes weighted by the phonon 
momentum (m). We note that the authors of Ref. [|J included a Coulomb term in this 
expression, which is incompatible with the lattice results. The flux tube model predicts the 
first excited T,' g + to be split by iV = 2 from the ground states, and two degenerate and 
E~ potentials at iV = 3. In the flux tube model the lowest II state is predicted to be the II U . 
It is split from the qq ground state by N = 1 and is followed by the N = 2, U g potential. 
The two lowest A states, A g and A u correspond to N = 2 and N = 3 respectively. 

The flux tube model fits the first excited state, U u quite well over a wide range of the 
quark separation. It is, however, the only surface to do so at small distance. Furthermore, 
this may be a fluke due to the particular choice of the short distance cutoff of the tt/R 
term employed in Ref. ||. Alternatively, at large distances, the system is expected to act 
like some sort of string with an excitation energy given by tt/R. The splittings do indeed 
appear reasonable for all the surfaces considered. It is, however, disconcerting that the II 
surfaces diverge from the flux tube model predictions for r > 4i? . This must be taken as 
an indication that the simple flux tube model considerations fail for more complex gluonic 
configurations. 

In Fig.3a-c we plot the results of our calculations for S, II and the A potentials with 
the flux tube results (solid lines) for comparison. It is apparent that the U g surface lies 
below the U u surface while the A u lies below the A g . This is opposite to the lattice results, 
indicating that the model fails to reproduce the expected level orderings. The £ levels are 
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also permuted with respect to the lattice calculations at small R. We note, however, that 
the correct level orderings may be reinstated by simply flipping the intrinsic parity of the 
gluon (set rjg — 1 in the expressions above). The resulting surfaces agree reasonably well 
with the lattice at small distances. For example, the U u and U g potentials are roughly 1 
and 2 at the origin in both calculations. We obtain values of approximately 3 and 3.5 for 
the A u and A g surfaces at the origin, similar to the lattice values of 2 and 3 respectively. 
The agreement persists to intermediate r, after which it is apparent that the model results 
approach the linear regime much too quickly with respect to the lattice. This is a strong 
indication that more degrees of freedom become active as the quark separation increases 
beyond r ~ R Q ~ 1/2 fm. This is certainly sensible from the flux tube model point of view, 
where a large number of degrees of freedom are necessary to construct string phonons. 

Our results include the contribution due to direct interaction between the static quarks 
(the last term of Eqn. (5)). This term was set to zero by Horn and Mandula |7j because 
it anti-confines and therefore may be expected to not produce a flux tube between the 
quarks. However, we note that it is responsible for producing a hybrid potential slope equal 
to that of the qq ground state potential, in keeping (roughly) with the lattice data. Thus 
we have chosen to retain this term. Note however, that this implies that there is a short 
distance repulsive Coulomb spike which should appear at very short quark separation. The 
appearance of this spike is, however, unphysical because the hybrid may emit a gluon and 
convert into a qq color singlet and a low lying glueball and this is energetically favorable for 
small r. Equivalently, Vu u > V q q + m gh for r < 0.2 fm. Thus the Coulomb spike mentioned 
above is irrelevant and one should see, instead, a Coulomb core for small r. This effect 
may be easily incorporated into the model presented here by allowing for coupling to the 
qqgg channel. Notice that no core is visible in the lattice data (especially in U u , which is 



measured down to r m 0.04 fm). The authors of Ref. j0| are currently examining this issue 
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Fig. 3b. II surfaces. Model (symbols) and flux Fig. 3c. A surfaces. Model (symbols) and flux 
tube (lines). tube (lines). 

III. CONCLUSIONS 



The model presented here agrees moderately well with lattice calculations of adiabatic 
potential surfaces. However, it disagrees in detail. In particular, it is necessary to ignore 
the intrinsic parity of the gluon to obtain the expected level orderings. Although one may 
argue that this is in keeping with the expectations for the parity of a lattice link operator 
(such as is employed in the development of the Flux Tube model), the lack of consistency 
is disconcerting. Furthermore, the model predicts surfaces which become linear too quickly 
with respect to the lattice. We expect this flaw to persist in all constituent glue potential 
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models (such as [0] and fL8f ). This is a strong indication that more degrees of freedom are 
required to describe soft glue at large interquark separation. 

It is important to note that models which employ single spinless gluons do not contain 
sufficient degrees of freedom to reproduce the potential spectrum. In particular, "single 
bead" Flux Tube models ( cannot make II g or E u states while three dimensional bead 
models ( 0) cannot make S + states. Thus including gluon spin is a minimal necessity in 
this class of models (although the level ordering problem must be overcome as discussed 
above). We also note that it is possible for spherical bag models to reproduce the lattice 
calculations at small quark separation, but that they fail at large R. Furthermore, Flux 
Tube model or Nambu-Goto string models [22| reproduce the lattice reasonably well for 



intermediate to large quark separations, but do not perform well at small R or in detail at 
large R. In particular the potential separation does not appear to be the expected tt/R and 
the slopes do not appear to agree. Finally, the bag model of Ref. @,^1|] works reasonably 



well over all quark separations, although problems remain to be resolved in the £ states and 
for large R. 

In summary, it appears that some sort of flux tube is required to explain the lattice 
adiabatic hybrid potential surfaces at large quark separation. This is in keeping with the 
conclusions of Ref. [^], where flux tubes were required to explain the spin splittings in heavy 
quarkonia. However, any flux tube model must attempt to incorporate the small R behaviour 
of the potential surfaces (and the seemingly anomalous behaviour at large R). This appears 
to be an indication that different degrees of freedom are required at small distances (as is 
incorporated in, for example, the bag model of Ref. ||). Finally, the model presented here 
(with the gluon parity reversal) works moderately well in describing the potential surfaces 
and could provide a useful starting point to simple models of gluonic hadron properties (note 
that the early transition to a linear potential seen in Figs. 3 should not be relevant since 
the hadron wavefunction is exponentially suppressed in this region). Thus we expect that 
the success of the previous glueball spectrum calculation (|TJ| in this model (where the gluon 
parity was irrelevant) was not a fluke. Several benefits of the model are particularly relevant 
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to hybrids; light quarks may be easily incorporated, the effects of coupled channels may be 
examined, and the effects of light quark coupling to virtual transverse or Coulomb gluons 
may be included. The latter effect is of interest because it is excluded in quenched lattice 
calculations and may have a substantial effect on the hidden flavor hybrid spectrum. 
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